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ABSTRACT 

We report the detection of CO molecular line emission in the z=4.5 millimeter- 
detected galaxy COSMOS_J100054+023436 (hereafter: JlOOO+0234) using the 
IRAM Plateau de Bure interferometer (PdBI) and NRAO's Very Large Array 
(VIA). The ^2C0(4-3) line as observed with PdBI has a full line width of ~ 
1000 kms~^, an integrated line flux of 0.66 Jy kms~^, and a CO luminosity of 
3.2x10^° Lq. Comparison to the 3.3a detection of the C0(2-l) line emission 
with the VLA suggests that the molecular gas is likely thermalized to the J=4-3 
transition level. The corresponding molecular gas mass is 2.6 x 10^° M0 assuming 
an ULIRG-like conversion factor. From the spatial offset of the red- and blue- 
shifted line peaks and the line width a dynamical mass of 1.1 x 10^^ Mq is 
estimated assuming a merging scenario. The molecular gas distribution coincides 
with the rest-frame optical and radio position of the object while being offset by 
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0.5" from the previously detected Lya emission. JlOOO+0234 exhibits very 
typical properties for lower redshift (z~2) sub-millimeter galaxies (SMGs) and 
thus is very likely one of the long sought after high redshift (z>4) objects of 
this population. The large CO (4-3) line width taken together with its highly 
disturbed rest-frame UV geometry suggest an ongoing major merger about a 
billion years after the Big Bang. Given its large star formation rate (SFR) of 
> 1000 M0?/r~^ and molecular gas content this object could be the precursor of 
a 'red- and- dead' elliptical observed at a redshift of z=2. 

Subject headings: galaxies: formation - galaxies: high-redshift - galaxies: ISM - 
galaxies: starburst - submillimeter 



1. Introduction 



Wide field blank sky surveys at millimeter and sub-millimeter wavel engths have estab - 
lished a population of active star forming galaxies at high redshift (e.g. iBlain et al.ll2002l ). 
These sources or so-called submillimeter galaxies (SMGs) dominate the (sub-)millimeter 
background, and represent (50-75)% of the star formation at high redshift caus ing a signif- 



icant revision to the opt i cally derived star f ormation history of the Universe (iSmail et al. 
I997I : iHughes et aPllQQsl : iBarger et allllQQsl : iBertoldi et aPboOoh . 



The bulk of this population lies at redshifts below z=3, with far- infrared (FIR) lumi- 
nosities of > 1O^^L0, which (if powered by star-formation) imply star-formation rates in 
excess of 1000 M© yr~^. This is sufficient to build up the stellar mass of a giant elliptical 
galaxy in about 1 Gyr (jChapman et al.ll2005l ). At redshifts z<3 about 100 SMGs are now 
known, many of which have been stu died in detai l , incl uding CO observations that indicate 
that they are very massive svstems JCreve et al. 2005[) . Rece nt high resolution studies of 



the molecular gas in z~2-3 SMGs (iTacconi et al.l 120061 . |2008| ) showed that the star form- 
ing regions are fairly compact and that the SMGs resemble "scaled-u p and more gas-rich 
versio ns" of the local Ultra-Luminous Infra-Red Galaxies (ULIRGs; e.g. iDownes &: Solomon 



19981 ). Due to their derived central densities, which are close to those of ellipticals or massive 



bulges, they appear to form stars at the maximal rate over very short time scales ("maxi 



murn starburst"). Thus the SMGs phase appears to last for about 100 Myr (ITacconi et al. 



20081 ). These massive starb urst galaxies place tight constraints on galaxy formation mod- 
els (e.g. iBaugh et al.ll2005l ). One major focus of current and future (sub-)mm surveys is 
to identify these massive starbursts at z>3, a redshift range for which SMGs could place 
tight constraints on hierarchical merger vs. monolithic collapse models in galaxy formation 
scenarios. 
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The discovery of such an SMG above a redshift of z=4 indicates that massive galaxy 
formation is aheady well under- way when the universe is only 1.5 Gyr old. The recently 
discovered o bject JlOOO+0234 originally selected as a V band drop-out with a weak radio 
counterpart (jCarilli et al.ll2008l ) and with a stellar mass of M^. > 10^° Mq, a young starburst 
age of 2-8 Myr and a star formation rate of SFR > 1 000 MtDjyr"^ is a c andidate for such a 
z>4 SMG, as it was detected in its FIR continuum (ICapak et al.l 120081 ). Throughout this 
paper we assume a standard concordance cosmology (Ho = 70, Qm = 0.3, = 0.7). 



2. Observations 

The bolometer camera MAMBO-2 at the IRAM 30m telescope was used in on-off mode 
in December 2007 and January 2008 to measure the 1.2 mm flux density of JlOOO-l-0234. The 
weather conditions during the observations were good, the reduction was performed using 
MOPSIC (written by R. Zylka, IRAM). The achieved rms for a total on-sky integration of 
70min is 0.67mJy. JlOOO-l-0234 was detected at a 5a level with a total 1.2mm continuum 
flux density of 3.4 mJy. 

The ^^CO(4-3) line tracing cold molecular gas in JlOOO-l-0234 was observed with the 
IRAM Plateau de Bure interferometer (PdBI) between January and April 2008 in the spe- 
cial C and C configurations. Both receivers were tuned to the redshifted line frequency of 
83.3257 GHz covering a total bandwidth of ~ 0.9 GHz. For calibration a nd mapping, we used 



the s tandard IRAM GILDAS software packages CLIC and MAPPING ( IGuilloteau fc Lucas 



20001 ). The quasar 1005-1-058 was used for atmospheric calibration while standard calibrators 
were used for fiux calibration. The 12 hr of total integration time result in a data cube of 
2.3" X 1.9" resolution (using robust weighting) with an rms of 0.42mJybeam~^ for 20 MHz 
(~ 72 kms~^) wide channels. 

We used the VLA in D configuration to observe the ^^C0(2-l) line in June and July 2008 
for a total of 21 hr. The two 50 MHz wide correlator units (each corresponding to a velocity 
width of ~ 360 kms'^) were tuned to 41.6352 GHz (on-line) and 41.7351 GHz (off-line). 
The observations were done in the fast-switching mode using the nearby quasar 1018-1-055, 
at a distance of 5.3° from our source, for atmospheric amplitude and phase calibration. 
20 antennas were available during the observations and the phase stability was typically 
~ 10°. The final images using natural weighting have a resolution of 1.9" xl.4" and an rms 
of 49 /iJy beam~^. 

All observed and derived properties of JlOOO-l-0234 are summarized in Tab. [H 



-4- 



3. The Molecular Gas Properties 



The CO(4-3) line was detected at a 5.5a level (Fig. [T]) in a 280 MHz wide channel, 
which corresponds to a full line width at zero intensity (FWZl) of ~ 1000 kms-\ i.e. with 
an integrated flux of 0.66 Jy kms~^. The peak position of the CO(4-3) emission (see Fig. 
[2]) agrees very well w 'ithin the pos i tional uncertainties with the rest-frame near-IR and radio 
position derived by ICapak et al.l (120081 ). The line emission of CO (4-3) is centered at an 
observed frequency of 83.1857 GHz, which corresponds to a redshift of z=4.5423 (assuming 
a rest frequency of 461.040768GHz for ^^CO(4-3)). Given the l arge line width th is is in 
agreement with the optically derived redshift of z=4.547±0.002 ( Capak et al.l [2OO8I ) . The 
CO (4-3) emission coincides with the rest-frame optical counterpart (as traced by the IRAC 
data) of JlOOO-l-0234 (Fig. [2]). Averaging the remaining channels blue- and redwards of the 
detected line emission results in a line-free image of 620 MHz width (~ 2200 kms~^) with 
a tentative 2.2a detection of the continuum at the rest-frame optical/radio position of the 
source (Fig. [T]). 

The VLA observations with a channel width of ~ 360 kms~^ are centered at the mid- 
dle of the line detected in CO(4-3). C0(2-l) line emission is detected at the 3.3cr level 
and coincides spatially with the CO (4-3) emission within its positional uncertainties of 
FWHMfeeam/(S/N)~0.5". Only about a third of the total line width was covered during 
the VLA observations, so the detected line flux of 0.057 Jy kms~^ could be higher by a 
factor of ~3, taking the missed blue- and redshifted emission into account by assuming a 
box-car line shape. The ratio of the CO(4-3) line flux of 0.66 Jy kms~^ to the thus corrected 
C0(2-l) line flux of ~0.16Jy kms~^ is about 4, implying that the molecular gas is still 
thermalized at the J=4-3 transition, as the flux still increases with z/^ at this transition. 

Continuum emission from J1000-f0234 was detected at 1.2 mm at (3.4±0.67) mJy, while 
the off-line channels in the interferometric observations resulted in a 3a upper limit of 150/iJy 
at 41.7GHz (VLA) and a tentative 2.2cr detection of ('0.17±0.0 8)mJy at 83 3 GHz (PdBI). 
All these values are consistent with the SED flts presented bv ICapak et al.l ( 2008o and a 



dust temperature of 30 - 50 K, typically found in z~2 SMGs (jPope et al.ll2006 : 
2006h . 



Kovacs et al. 



^The scaling of the SEDs presented in Fig. 4 of ICapak et al. 
(P. Capak, priv. comm.). 



pOOSf ) is overestimated by a factor of 10^ 
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4. Masses and Typical SMG Properties 



As the molecular gas appears to be thermalized at least up to the CO (4-3) tran- 
sition as inferred above, its line flux can be used to estimate the CO luminosity [set- 
ting Lqq(J = 4 — 3) = L'(^o(J = 1 ~ 0)] hence the molecular gas mass present in 
JlOOO+0234. To der ive the CO luminosity L^ q and molecular gas mass Mgas we use equa- 
tions 3 and 4 from ISolomon fc Vanden BoutI (120051 ) with a conversion factor for L^ n to 



-lpc2 



^ derived for local ULIRGs (jPownes fc Solomonlll998l ) and 



Mgas of a=0.8MQ(K kms' 
used for high redshift objects (jSolomon fc Vanden BoutI l2005l ). The observed line flux of 
Sco(4-3)=0-66 Jy kms"^ corresponds to a CO luminosity of L'(;.q ~ 3.2 x 10^° Lq or a molec- 
ular gas mass of Mgas ~ 2.6 x 10^° Mq. 



We imaged the red and blue (~ 500 km s~^ wide) halves of the line emission separately, 
resulting in a 4.7(j (red) and 3.8cr (blue) emission peak, respectively. The two peaks (Fig. 
E]) show a spatial shift of ~ 1" roughly from northwest to southeast (with a positional 
uncertainty of ~0.5"). This offset corresponds to ~ 6.6 kpc at the redshift of JlOOO+0234, 
suggesting that the CO emission might be fairly extended. Using this spatial offset between 
the red- and blue-shifted half of the line emission and the line width of ~1000 kms^^, we can 
estimate the dynamical mass of JlOOO+0234. We use the relation between dynamical mass 
Mciyn, velocity width Avfwhm = Vrot(r) x sin z/2.4 and ra dial extent r of the CO emitting 
region: Mdyn x sin^z = 4 x 10^ r A Vp^jjjyj ( iNeri et al.l 120031 ). As the system is likely merging 
(as discussed below) we include also a factor of 2 for a merging s ystem assuming that the gas 
has already settled into the new potential (iGenzel et al.ll2003l ). The estimated dynamical 
mass of JlOOO+0234 is about 1x10^^ Mq for values of Vrot(r) x sin z ~ 1000 kms~^ and a 
radius r of 3.3 kpc. Thus JlOOO+0234 has a gas fraction of Mgas/M^yn ~ 25%, typical for 
SMGs (iGreve et al.ll2005l : iTacconi et al.ll2006l ). We would like to caution that the derived 
offset of ~1" has signiflcant uncertainties and thus the derived dynamical mass could be 
overestimated if the separation were smaller. 

The properties of JlOOO+0234 are very similar to the global properties derived for z~2 
SMGs in terms of molecular gas mass (< Mgas >~ 3 x 10^° M©), extent (rco ~ 2 kpc), CO 
line width (< FWHMco >~ 780 kms"^), dynamical mass ( < Mdyn >~ 1-2 x 10" M©) 
and molecular gas mass fraction (~ 25%) (IGreve et al.ll2005l ). Its estimated current gas 



consumption tim e scale of t = Mga s/SFR ~ 30Myr assuming our gas mass and a SFR of 
> 1000 MQyr~^ (ICapak et al.ll2008) is si milar to the gas depletion rates of ~ 100 Myr found 
for the z~2 SMGs (E^coni et al.l |2008l). Taking together the stellar mass of > 10^° M© 



produced in a recent burst ~ 7 Myr ago (ICapak et al.ll2008l ) and the available gas reservoir 
and gas consumption time scale derived here, a signiflcant fraction of the stellar mass of a 
massive elliptical could be produced on a relatively short time scale. Assuming that these 
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recently formed stars will be passively evolving since then, JlOOO+0234 could turn into a 
'red-and-dead' elliptical - like those that are found at z=2 when about 2 Gyr have passed 
since z=4.5. 



5. Evidence for an Ongoing Major Merger 



The molecular gas emission arises from a region that is highly obscured at the observed 
optical wavelengths (see Fig. 1 of ICapak et al.ll2008l . Fig. [2] and [3]). However, it is coincident 
with the position of th e rest-frame optica l and radio emission (Fig. 2j). The HST ACS 
image of JlOOO-l-0234 JCapak et all boosh was adaptively smoothed dScoville et all boooh 
to enhance the structure present in the blue continuum covered by the F814W filter (Fig. 
[3]). The molecular gas appears to lie next to bright blue continuum emission that is at its 
eastern side. Lya e mission has also been associated with this bright rim in blue continuum 
(jCapak et al.ll2008l ). The motion of the gas is roughly along a region of higher extinction 
(compared to the blue western component) running roughly from north to south and shows 
a larger spatial extent than the blue continuum emission. Note th e object west of th e gas 
emission is a foreground galaxy at z=1.41. As already mentioned by lCapak et al.l (120081 ) this 
geometry is very reminiscent of an ongoing merger. As the redshift of the molecular gas is 
very close to the redshift derived from the Lya line (z=4.542 vs. z=4.547), it is very likely 
that both components belong to the same (highly disturbed) system. Both the large FWZI 
of the CO (4-3) line and the appearance of JlOOO+0234 in Lya make an ongoing merger 
scenario very likely. We derive a conservative estimate of the merging time for JlOOO-l-0234 
of significantly less than a billion years using the prescription of iKitzbichler fc Whitd (120081 ) 
for close galaxy pairs with velocities of f < 3000 kms^^. This time scale is consistent with 
having a relaxed system at a redshift of z=2. 



The expected number of z=4-5 SMGs from the model of iBaugh et al.l (120051 ) yields of the 
order of 10-20 such sources within the 2 deg^ COSMOS field . As JlOOO-l-0234 was originally 



identified as a Lyman Break Galaxy via its V-band drop (ICapak et al.l l2008l : ICarilli et al. 



20081 ) . this suggests that UV emission can be detectable from ongoing major mergers at these 



redshifts. We have identified 4 additional V-band dropout LBGs with 4cr radio detections 
(jCarilli et al.ll2008l ) that could also be SMGs at this redshift. An inspection of their optical 
appearance in the HST ACS data shows that two of them also have a disturbed geometry. 
This suggests that there could be more objects sharing properties of both the SMG and LBG 
population at z>4 and JlOOO+0234 might therefore provide an important link between the 
SMG and LBG populations at high redshifts. 



The case of JlOOO+0234 shows that (at least some) z>4 SMGs can have very similar 
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properties to their well-studied low redshift (z~2) counterparts, that major mergers with 
very large SFRs are likely present at z=4.5 and that these systems might be the precursors 
of elliptical galaxies found at the red sequence at a redshift of 
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Table 1. Properties of JlOOO+0234 



Property 


Value 


Comment 


S250GHZ [mJy] 


3.4±0.67 


MAMBO 


SssGHz [mJy] 


0.2±0.09 


off-line channels (~2(T, PdBI) 


S42GHZ [mJy] 


<0.15 


off-line channel (<3cr, VLA) 


ZC0{4-3) 


4.542 


center of CO (4-3) line 


FWZIco{4-3) [kms'i] 


~ 1000 


CO (4-3) line 


Sco(4-3) [Jy kms"^] 


0.66±0.12 


5.5cr detection (PdBI) 


Sco(2-i) [Jy kms"^] 


0.057±0.017t 


3.3cr detection* (VLA) 


R.A.co{4-3) (2000) 


10:00:54.484 


line peak 


Dec.co(4~3) (2000) 


+02:34:35.73 


line peak 




3.2x10^° 


from CO (4-3) 


Mgas[M0] 


2.6x10^° 


from CO (4-3) 


Mdyn X sin2^[M0] 


1.3x10" 


using FWZI & r=0.5" 



Note. — Some properties measured f rom our PdBI and VLA data and 
deri ved using t he eq uations given by ISolomon fc Vanden BoutI (120051 ) 
and lNeri et all J2OO3L 

^ observed line flux not corrected for small bandwidth (see text for de- 
tails) 

* available bandpass corresponds to Av?^360 kms~^, i.e. a third of the 
total line width 
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54.8^ lO" 0""54.5^ 54.2^ 54.8^ lo" 0""54.5^ 54.2^ 



R.A. (J2000) R.A. (J2000) 



Fig. 1. — CO(4-3) line emission and corresponding continuum in JlOOO+0234 as observed 
with the PdBI, shown on the same color scale, but with different contour levels. Continuum 
emission from JlOOO+0234 is detected at the 2a level, while the emission from the CO(4-3) 
line shows a solid detection at 5.5cr. The integrated CO(4-3) line emission {lejt) is centered at 
a velocity of +500 kms~^ relative to the observed frequency of 83.3257 GHz. Contours start 
at 2(7 in steps of \a with \o = 0.12 Jybeam"^ kms~^. The integrated continuum (right) has 
the same contour steps, however la corresponds to 0.08 mJy/beam. The cross indicates the 
position of the CO (4-3) peak. The beam is shown in the lower left corner. 
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R.A. (J2000) R.A. (J2000) 



Fig. 2. — Location of the molecular gas relative to other tracers. Left: The CO (4-3) line 
intensity overlaid in contours (same as in Fig. [1]) on the Subaru ip band image, showing an 
offset from the bright blue continu um rim. Right: VLA 1 .4 GHz continuum contours (65% 
to 95% of the maximum of 45/iJy; ISchinnerer et al.l 120071 ) overlaid onto the CO(4-3) map 
showing an excellent correspondence between the molecular gas and the radio continuum, 
which is presumably tracing star formation activity. 
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Fig. 3. — Integrated red- and blueshifted CO (4-3) line emission in contours overlaid on the 
adaptively smoothed HST ACS/F814W image of JlOOO-l-0234. A motion from southeast to 
northwest is apparent with an offset of ~1" between the peaks of the red- and blueshifted 
emission. The contours start at 2(j in steps of la with la = 0.085 Jybeam"^ kms~^ and 
0.084 Jybeam"^ kms~^ for the red- and blueshifted line emission, respectively. The object 
~ 0.5" west of the CO emission is a foreground galaxy at z=1.41 (jCapak et al.ll2008l ). 



